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ABSTRACT 

Phase distortion due to compressible, inviscld flow over small perturbation 
laser turrets in subsonic or supersonic flow was calculated. The turret shape was 
determined by a two-dimensional Fourier series; in a similar manner, the flow 
properties are given by a Fourier series. 

Phase distortion was calculated for propagation at several combinations of 
elevation and azimuth angles. A sum was formed from the set of values, and this 
sum became the objective function for an optimization computer program. The 
shape of the turret was varied to provide minimum phase distortion. 


INTRODUCTION 

For many applications of a high energy laser on board an aircraft, the beam 
must be propagated with minimum phase distortion. The well known Strehle 
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relation [1] gives the decrease in far field intensity as a consequence of the rms 
phase distortion. As a result of the dependence of index of refraction on mass 
density, compressible flow over laser turrets causes phase distortion. The variable 
density and variable index of refraction surrounding an aircraft laser turret may 
be the result of viscous flow phenomena or inviscid flow. This paper focuses on the 
inviscid flow problem. 

In regard to the solution of the phase distortion problem due to inviscid 
compressible flow, several options exist. Location of the turret on the aircraft 
is an important consideration. Adaptive optics may be used. Correct shape of 
the turret can reduce significantly phase distortion. The approach taken in this 
paper is to minimize phase distortion through turret shape. 

An alternate approach would be to consider combined adaptive optics and turret 
shape. Higher order distortions, e.g., astigmatism and coma, are more difficult 
to remove by adaptive optics than lower order distortions. Using this approach to 
design, the turret shape should be modified so as to minimize higher order 
distortions. 

In passing, a comment should be made about adaptive optics. Adaptive optics 
for compensation of atmospheric turbulence and thermal blooming requires 
mirror displacements of a fraction of wavelength at frequencies of 25 kHz or so. 

In contrast, the adaptive optics for compensation of laser turret phase distortion 
requires mirror displacements of a few wavelengths at frequencies of a few Hertz. 

An analytical model for describing laser turret geometry and the associated 
compressible flow field has been described in the papers by Fuhs [2] and Fuhs and 
Fuhs [3,4]. A companion paper in this conference proceedings [5] discusses the 
analytical model. The turret shape is described by a two-dimensional Fourier series. 
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Using the flow over a wavy wall on a circular cylinder as the basic solution, 
a Fourier series can be found for the potential function. In contrast to a 
direct numerical integration of the equations of motion of gas dynamics, an 
explicit analytical solution is obtained. As a result, computer time is significantly 
less. 

At a plane normal to the beam sufficiently far from the aircraft, the phase 
distortion is calculated. One method to represent the phase distortion is to use a 
series with Zemike polynomials [6]. The advantage of using a Zernike series is 
that the coefficients in the series are related directly to the magnitude of the 
various types of distortion, i.e., tilt, focus, defocus, astigmatism, etc. 

A computer program [7,8] has been developed which can find the values of design 
variables yielding a minimum value of an objective function subject to constraints. 

The program has been applied to a variety of problems [9,10]. For the case at hand, 
the design variables are the coefficients describing the turret geometry. Also 
gas density inside the turret and the location of the laser turret primary mirror 
were treated as design variables; and in Figure 1 define location of the 
mirror. The objective function was the weighted sum of phase distortion for several 
sets of elevation and azimuth angles. Constraints included the maximum slope of the 
laser turret as well as maximum turret height. 

PHASE DISTORTION 

Optical path length, L^, is defined as 

b 

= j* n(s)ds (1) 

a 

where n is the index of refraction and s is distance along a particular ray. Points 
a and b are on the ray. The subscript i identifies the ray. The difference in 
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Liepmann and Roshko [11] as 
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( 5 ) 


C = - 2u - 
P 

where u is the perturbation velocity in the freestreain direction and v is the 
radial velocity which is normal to the fuselage axis. Equations (1) to (5) can be 
combined to give an integral for the phase distortion in terms of pressure coefficient; 
see equation (7) of Reference [5]. When the potential function for the flow is known, 
both u and v can be calculated. 

METHODS TO COMPENSATE FOR PHASE DISTORTION 

The options available for compensating for phase distortion due to inviscid 
compressible flow over the laser turret were mentioned in the Introduction. Adaptive 
optics is one technique. Wolters and Laffay [12] demonstrate the effectiveness of 
adaptive optics. 

Another method to compensate for phase distortion is to use a laser turret 
of proper shape. This is the approach of this paper. Turret geometry constitutes 
a passive technique. Turret geometry as a means to lessen phase distortion 
is discussed in the following sections. 

A method to represent phase distortion is to use Zernike polynomials [1,6]. 

The polynomials, which are given in the paper by Hogge and Butts [6], are an 
orthonormal set of functions. The phase distortion is 

n 

P = Z AjFj (6) 

j=l 

where A, is a coefficient and F. is the jth Zernike polynomial. The summation 
extends from 1 to n, where F^ is the highest order polynomial considered. Typically 
n = 10 is adequate. The coefficient A^ is obtained by multiplying both sides of 
equation (6) by F and integrating over the aperture or beam cross section. All 

q 
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terms in the summation vanish except for the term j = q. Equation (19) of 

Reference [5] gives the results and a formula for A.. 

3 

An alternate method to compensate for phase distortion is to combine adaptive 
optics and turret geometry. Higher order phase distortions are more difficult 
to compensate by adaptive optics. Form an objective function which is 


B = y W.A. 
^ 1 1 


(7) 


where is a weighting factor for the ith coefficient in Zernlke series for phase 
distortion. The larger i, the larger is the value of W^. As the turret shape is 
varied, the value for B changes. Using COPES/CONMIN computer program [7], the 
value of B can be minimized through variations of turret geometry. The consequence 
is that the effectiveness of adaptive optics is enhanced since higher order 
distortions are minimized. 

The technique of adaptive optics employs segmented or deformable mirrors. 

Phase is controlled by mirror displacement. Further, the frequency of the mirror 
motion is determined by the frequency of the adverse phenomenon being overcome 
through use of adaptive optics. Hence, different types of adaptive optics can be 
thought of as occupying different locations in the mirror amplitude/f requency plane. 
Compensation for atmospheric turbulence occurs in the low amplitude, high frequency 
region of the adaptive optics map. Compensation for the adverse influence of flow 
over the turret occurs in the high amplitude, low frequency region of the adaptive 
optics map. The frequency response is dictated by turret slew rates or aircraft 
maneuver rates. 
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The (+) sign is for subsonic flow, and the (-) sign is for supersonic flow. 

The quantity 3 is 

6" = |1 - M^l (9) 

where is the freestream Mach number. The potential (f> is the perturbation 
potential and yields the perturbation velocities 


3x ’ 3r ’ r 99 


( 10 ) 


The velocities appearing in equation (5) can be obtained from equation (10) . 

The boundary conditions for equation (8) and a wavy wall on a circular cylinder 
are discussed in Reference [4]. A solution is obtained for one spatial frequency 
for the wavy wall. The solution is the basic function from which a Fourier series 
for the flow is constructed. The turret is represented by two polynomials 

K 

Z 'Xf k 


( 11 ) 


k=l 


and 


P 

f(6) = 1 + y 

j 


( 12 ) 


To obtain a symmetric turret in the Q-direction, only even values of j are used in 
equation (12). In terms of f(x) and f (0) , the turret geometry is 

R(x,0) = Ro + ef(x)f(9) (13) 

where R is the radial distance to the surface of the turret or fuselage. 

Equation (13) is represented by a Fourier series which leads to the coefficients 
in the Fourier series for the potential flow. 

Figure 2 shows the geometry. The maximum turret height is e, and the length 
of the turret is 2K'. The meaning of 9 is that f(0) is zero for |0l > 6 
Figure 3 is an artist’s concept of the laser turret. Two comments are applicable to 
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Figure 2. Geometry of a Small Perturbation Laser Turret on a Circvilar Fuselage. 



Figure 3. Artist's Concept of Laser Turret. 

Figure 3. First, the turret will become operational at some future date when laser 
canopies can be manufactured I Second, the model does not include perturbation effects 
of wings, blade antennas, and similar items. 
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Figure 4 illustrates the coordinate system used to describe the direction of 
the beam relative to the aircraft. A Cartesian coordinate system X, Y, Z is 
oriented as shown. The Z-axis forms the polar axis for a spherical coordinate 
system. The beam is at azimuth angle (}) and elevation y* 



Figure 4. Coordinate System for Direction of Laser Beam Propagation. 


II II 

The design variables become a, and b.. Conditions at |9[ =0 and lx[ = £ 

iC J IQ^X 

reduce the number of independent design variables. Typical conditions are 

(14) 


R(6,£) = Ro 

max 


and 


9R) 

8xj 


= t) 

■/x=+£ 


= 0 


(15) 
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The intercept between a ray and the turret surface is a particularly difficult 
problem in analytical geometry. An iteration scheme was used to find the intercept 
as follows : 

a. At a given point s = s’ on the ray, calculate x,0,R. 

b. For these values of x and 0, calculate the radial component of the 
turret surface R^. 

c. Iterate to find a value of s' such that R = R^. 

The iteration procedure is described in detail by Vanderplaats [13]. 

In equation (1), the integration starts at a and ends at b. The start of 
integration, a, is at the primary mirror surface. The end of integration, b, is at a 
point sufficiently far from the laser turret so that additional integration by an 
amount 6s yields a negligible change in the value of the integral. Since integration 
starts at the mirror surface, a portion of the ray between a and b is within the 
laser turret. Hence the density within the turret is a factor in the phase 
distortion. Even if the air external to the turret were uniform, a phase distortion 
could be generated by the air within the turret. In this study, the turret 
window is assumed to be of zero thickness so as to be distortionless. 

OPTIMIZATION OF LASER TURRET SHAPE 

For any given azimuth, (j), and elevation angle, Yj the phase distortion of any 
ray in the beam can be calculated using the center ray of the beam as a reference. 

At a specified beam orientation, the phase distortion typically will be calculated 
using equation (7) of Reference [5] at two radial locations for each of eight 
angular locations, i.e., r) of Figure 1 occurs every 45°. Furthermore, to provide 
optimum overall system performance, several beam orientations are considered. By 
squaring the values of phase distortion and summing over all rays and orientations, 
a measure of total performance, S, is obtained as 
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= 2 Z V 2 Z 

(Pr z T) 


(16) 


The variable z is defined in Figure 1; z gives the radial location within the beam. 

is a weighting function. Values of are determined from mission studies. 

For a particular mission, the laser beam may be pointed most of the time at a 


())Y 


particular direction, i.e., particular values of (p and Y- For that direction, W 
is larger. For more extensive mission studies, sufficient information can be 
obtained so that a meaningful function can be defined for the function is 

a two-dimensional probability density function giving the probability the beam 
points in the direction specified by <f) and y. 

The objective of the optimization was to minimize S by determining the 
proper combination of design variables. The design variables have been 
mentioned earlier. To summarize, the design variables are as follows: 

independent variables from 
independent variables from 
mirror location 
density within turret 

The analysis capability presented in this paper and Reference [5] has been 
coded in FORTRAN to produce maps of phase distortion for given azimuth and 
elevation angles. Figures 5 to 7 are phase distortion maps. For all three maps 
the windward side of the beam is at the top of the map. This fact can be determined 
from the equation for phase distortion 


O/ 
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j (Sj - s^)p^ J 


yM^C ^ 


_E] is _ 

2 'j Rq Po 


I 


yM^C ^ , 
--fLJB-i is 
2 JiR^ 


(17) 


s , 
3 


s . 

X 


349 




In equation (17) , Is the distance from the surface of the primary mirror to the 
intercept of ray j with the laser turret surface. The laser radiation has wavelength 
X. When the pressure coefficient along ray j exceeds that along ray 1, a positive 
contribution to the value of phase distortion occurs. Refer to Figure 2 or Figure 4. 
For an elevation angle of y = 45° and the laser beam In the plane of symmetry of the 
fuselage, s > s . Consequently, the term (s. - s.)P^ in equation (17) is also 
positive on the windward side. A positive phase distortion, P, means a lag of 
wavefront for ray j compared to wavefront of ray 1. 

The laser turret geometry and the associated flow were coded in subroutine form 
compatible with the general purpose optimization program. 


COPES /CONMIN [7] 

The COPES/CONMIN program solves the design problem of the following form: 

Minimize F(x) (18) 

subject to the constraints 


Gj (x) ^0 1 j = 1 . . . .m 


(19) 


where F(x) is called the objective function and is defined by equation (16). 

The vector of design variables, x, contains the design variables summarized above. 

G.(x) are the constraints. The constraints that were considered at one time or 
3 

another during the study were the maximum slope of the turret in the streamwise 
direction, no discontinuity at turret fuselage boundary, and the conditions of 
equation (15). The slope was restricted, at the most, to a value of 0.3 since 
the linearized flow equations become inaccurate for larger values. The choice of 
objective and constraint functions is somewhat arbitrary; the only restriction is 
that both functions must be continuous functions of the design variables, x, with 
continuous first derivatives. In general, F(x) and G^ (x) may be any linear or non- 
linear functions of x. 
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TWO EXAMPLES OF LASER TURRET DESIGN 


Two design examples are presented here, the first being for subsonic flow 
and the second being for supersonic flow. The design conditions are listed in 
Table I. Calculations were conducted for six beam orientations and sixteen rays 
within the beam. 


Table I. Design Conditions 


AERO-OPTICS 


Mach number Case 

1 

M^ = 0.5 

Case 

2 

M =2.0 
00 

Ratio of heat capacities 


Y = 1.4 

Wavelength of laser radiation 


A = 3.8 microns 

Density ratio 


P=./PSL ' 

Constant for index of refraction 


<’ = 0.00023 


GEOMETRY 


Fuselage radius 
Spacing of turrets 
Mirror location 

Turret length 
Turret height 
Mirror radius 

Maximum angle extent of turret 


L = 5.0 
"m = "-125 

I = 2.0 


E = 0.2 

R =0.05 
m 

e = 60° 

max 


GEOMETRIC BOUNDARY CONDITIONS 


x / % 

f(x) 

df (x) / dx 

0/e 

max 

f(6) 

df(0)/d9 

1.0 

0 

0 

+ 1.0 

0 

varies 

0 

1.0 

varies 

0 

1.0 

0 

1.0 

0 

0 
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Table I Continued. Design Conditions 



BEAM ORIENTATIONS 


Number 

Azimuth, (f), degrees 

Elevation, y, degrees 

1 

0 

45 

2 

0 

90 

3 

0 

120 

4 

45 

45 

5 

90 

30 

6 

90 

60 


PHASE DISTORTION CALCULATION POINTS 


Rays defined by all combinations of: 

radius within beam ^^^0 ~ 0*025, 0.050 

angle within beam U = 0, 45, 90, . . . 315 

Note: All rays are shown in Figure 1. 


CONSTRAINT IN SLOPE 


- 0.3 < < 0.30 

— dx — 

for 0=0 


As a reference, flow over a cosine-shaped turret was calculated. The equations 
for turret geometry were 


2 4 

f(x) = 1.0 - 0.50(~) + 0.0625(j) 


and 


f(0) = 1.0 - 1.824(|- ) + 0.832(| ) 

max max 


( 20 ) 

( 21 ) 
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Phase distortion maps are shown in Figures 5 and 6 for the turret specified by 
equations (20) arid (21). The perturbation velocities, u and v, were calculated 
in the plane of symmetry of the fuselage and are shown in Figures 8 and 9. Figure 8 
is for the subsonic flow example, and Figure 9 is for the supersonic flow example. 
The radial perturbation velocity, v, is dictated by the boundary condition 
at the turret surface. Hence v is identical in both Figures 8 and 9. The axial 
perturbation velocity, u, is different for subsonic flow as compared to supersonic 



Figure 8. Perturbation Velocities for 
a Cosine Turret in Subsonic 
Flow. 


Figure 9. Perturbation Velocity for 
a Cosine Turret in Super- 
sonic Flow. 


flow. In subsonic flow, the maximum value of u occurs at x = 0, while in 
supersonic flow, at x = 0, the value of u is zero. In supersonic flow, the flow 
is compressed (u < 0) on the forward or windward side of the turret; on the 
leeward side of the turret, the flow is expanded (u > 0) . 
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Results for the two examples are summarized In Table II. In Table II the values 

'll 'll 

of the coefficients a and b are listed for both the initial and optimized laser 

f\j 

turrets. Using these values of a and b, the laser turrets have been drawn and 
appear In Figures 10 to 12. Figure 10 is the cosine-shaped turret used as reference. 


Table II. Summary of Laser Turret Design Results 



Quantity 

Initial 

Turret 

Optimized 

Turrets 


Subsonic 

Supersonic 

Case 1 (Subsonic) 

Case 2 (Supersonic) 

S 

36.02 

2.69 

31.22 

1.55 

Pt^Pco 

0.7 

0.3 

0.7 

0.3094 


355 











Figure 12. Optimized Laser Turret for Supersonic Flow. = 2.0. 

Figure 11 is the optimized laser turret for subsonic flow. Figure 12 is the optimized 
laser turret for supersonic flow. Comparing Figures 10 and 12, very little difference 
between the cosine-shaped and the optimized turret for supersonic flow can be 
seen. However, reference to Table II shows the optimized turret has odd powers 
for f(x); note that = 0.2651, a^ = - 0.1326, and a^ = 0.0166. 

For the calculations summarized in Table II, all the weighting values 
were unity. Using the COPES/CONMIN optimization computer program, the objective 
function, S, defined by equation (16), was reduced from 36.02 to 31.22 for the 
laser turret optimized for subsonic flow. The reduction is 13 per cent. The objective 
function, S, for the laser turret designed for supersonic flow was reduced from 
2.69 to 1.55. The reduction is 42 per cent. Also note that a density within the 
turret, p^, less than ambient helps to reduce the phase distortion; see Table II. 
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EXAMPLE OF LASER TURRET GIVING LEAST AND WORST PHASE DISTORTION 


To illustrate the range of values of the objective functions, S, that can be 
obtained by varying the turret geometry, both the best and worst laser turret were 
designed. Details are given in Table III. The cross sections of the initial turret, 
the best turret,' and the worst turret are shown in Figure 13. For the case at hand, 
the initial turret yielded S = 0.0115. The worst turret gave S = 0.0918 which is a 
change of 690 per cent. The best turret has S »= 0.0012 which is an improvement in S 
of 890 per cent. The range from the worst to the best is 0.0918/0.0012 *= 78.5. 

Table III has the coefficients A^ for the Zernike polynomials. The phase 
distortion can be represented by equation (6) using A^ from Table III. The 
reader should compare A^ for the Initial turret with the other two turrets. The 
best turret has a slightly larger value for A^ . The value of Ag is reduced greatly. 

For the laser turret giving the worst distortion, all coefficients are increased 
except for A^. The focus coefficient is slightly smaller. The average value A^ 
has little significance. 

Compared to the two examples of the previous section, the turrets in this 
section were optimized for only one beam direction. The beam direction was at an 
azimuth, 4> = 45^ and an elevation, y = 45°. 

COMPUTER CODE FOR LASER TURRET OPTIMIZATION 

An extensive and versatile computer code has been written by Vanderplaats and 
Fuhs [14]. The computer code is based on References [2,3,5,7,8,13]. The program 
calculates the optical path length and phase distortion arising from the density 
field surrounding a laser turret in compressible flow. Further, the program finds 
the optimum turret shape yielding minimum phase distortion. The optimization 
and control codes are thoroughly discussed in Reference [14]. Sample data input and 
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Table III. Summary of Results for Phase Distortion 


Zernike 

Coefficient 

Physical 

Significance 

Initial 

Turret 

Laser Turret Shape 
Least Worst 

Distortion Distortion 


average value 

-3.331E-04 

3.967E-04 

-4.193E-03 

*2 

x-tilt 

0.02129 

4.391E-03 

0.05995 


y-tilt 

-7.643E-03 

5.422E-03 

>■0.03595 


focus 

-3.947E«04 

4.759E-04 

-3.650E-03 

^5 

astigmatism 

6.576E"05 

3.320E-05 

3.748E-04 

A 

6 

astigmatism 

9.183E-04 

3.973E-05 

9.944E-04 


coma 

-1.638E-05 

-7.172E-06 

-5.610E-04 


coma 

-5.123E-06 

3.136E-05 

-8.846E-04 


coma 

0.01656 

3.428E-03 

0.04549 

^10 

coma 

-5.951E-03 

4.285E-03 

-0.02840 


Mach niimber = .500 
Flight altitude = sea level 
Turret height /fuselage radius = .200 
Beam radius /fuselage radius = .05 
Elevation angle = 45° 

Azimuth angle = 45° 


sample output are given. The material is presented in sufficient detail so that 
Reference [14] constitutes a user’s manual for LASTOP. 
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axis of fuselage 


(a) Nominal turret. Objective function has value of 0.0115. 



(b) Turret yielding greatest distortion. Objective function has value of 0.0918. 



(c) Turret yielding least distortion. Objective function has value of 0.0012. 


Figure 13. Cross Section Shape of Turrets in the Plane of Symmetry. 
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COMMENTS AND CONCLUSIONS 


A versatile analysis and computer program has been developed which optimizes 
laser turret geometry to obtain minimum phase distortion. The turrets are located 
on a fuselage of circular cross section. Turret slope is limited so as not to 
exceed the perturbation allowed by the linearized equations for the flow. The 
computer code is described in Reference [14]. 

Examples have been given which show the decrease in objective function, S 
in equation (16) , that can be achieved. 
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